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I.   INTRODUCTION 

One  of  the  important  problems  in  the  theory  of  graphs,  especially  in 
its  applications  to  physical  systems,  is  that  of  locating  the  simple  cycles 
in  a  directed  graph.   The  author  conducted  an  extensive  literature  search 
(all  American  Mathematical  Society  Reviews,  I0E0E.E.  Circuit  Theory  Re- 
views, and  American  Computing  Machinery  Reviews),  and  was  unable  to  find 
any  really  efficient  methods  for  solving  the  problem.   One  of  the  more 
efficient  methods  was  the  algorithm  of  Danielson  [3],  which  uses  matrix 
techniques  and  is  designed  for  application  to  undirected  graphs,  although 
it  lends  itself  to  directed  graphs  if  certain  modifications  are  made.   In 
comparison  with  the  algorithm  which  we  shall  develope,  it  will  be  seen 
that  it  is  itself  rather  inefficient. 

We  will  first  give  a  number  of  fundamental  concepts,  with  examples, 
in  order  to  familiarize  the  reader  with  the  basics  involved,  and  in  order 
to  make  the  paper  as  nearly  self  contained  as  possible.   In  this  regard, 
we  have  adapted  considerable  material  from  Busacker  and  Saaty  [2]. 

After  a  consideration  of  the  basic  theory,  we  develope  an  algorithm 
for  locating  the  simple  cycles  in  a  directed  graph.   In  doing  so,  we  state 
some  specific  definitions  which  may  not  be  found  in  the  literature,  but 
which  are  fundamental  to  our  algorithm.  After  presenting  the  algorithm, 
we  give  an  example  to  show  how  it  works,  and  then  we  prove  that  the  algori- 
thm does  indeed  do  what  it  is  supposed  to  do. 

Following  the  statement  of  the  algorithm,  we  give  modifications  which 
allow  the  algorithm  to  be  applied  to  coalesced  graphs  (graphs  in  which  a 
subgraph  is  represented  by  a  single  vertex.   We  then  give  an  example, 
which  is  taken  from  circuit  theory,  and  show  how  the  algorithm  is  applied. 


II.   FUNDAMENTAL  CONCEPTS 

Let  S  be  a  set.   Suppose  s  is  an  element  of  S,  and  that  t  is  an  ele- 
ment of  S.   By  the  symbol  (s&t)  we  mean  the  unordered  pair  s  and  t  (un- 
ordered in  the  sense  that  (s&t)  and  (t&s)  denote  the  same  element).   We 
shall  use  the  symbol  (S&S)  to  denote  the  unordered  prod'  ct  of  S  with  it- 
self.  Now,  let  V  be  any  set  which  is  non-empty,  and  let  E  be  a  set  dis- 
joint from  V.   Define  a  mapping  g  from  E  into  (V&V) .   An  undirected  graph 
is  the  triple  (V,E,g).   We  denote  such  a  graph  by 

G  =  G(V,E,g), 
or,  simply 

G  =  G(V,E) 

if  the  mapping  is  understood.  We  call  the  elements  of  V  vertices,  the 

elements  of  E  edges,  and  g  is  called  the  incidence  mapping  for  G. 

A  directed  graph  is  the  triple 

G  =  G(V,E,g'), 
where  V  and  E  are  defined  as  above  and  where  g'  maps  E  into  (VxV) ,  the 
unusual  Cartesian  product  of  V  with  itself. 

We  use  the  terminology  of  Harary  [4]  and  call  a  directed  graph  a  di- 
graph. As  in  the  undirected  case,  we  call  the  elements  of  V  vertices  and 
the  elements  of  E  edges.   We  shall  denote  both  elements  of  V  and  elements 
of  E  by  subscripted  lower  case  letters.   If  g'  is  one-to-one,  we  shall 
sometimes  denote  an  element  e  of  E  by  the  pair  (v  ,v0)of  vertices  for  which 
g'(e)=  (v-,v0)-   Note  that  to  each  digraph  there  corresponds  an  undirected 
graph  formed  from  the  digraph  by  removing  the  direction  from  all  edges. 
A  ;  aph 

H  =  H(VH,ER,gH) 


is  said  to  be  a  subgraph  of  a  graph  G  =  G(V,E,g)  if  and  only  if  the  following 
conditions  hold: 

(i)  vH9  V, 

(ii)  EH£  E, 
(iii)  gH  =  g|H  , 
where  g   ,  denotes  the  restriction  of  g  to  the  subset  E„  of  E. 
If  e  is  an  edge  in  a  graph  G  =  G(V,E,g),  and  if 

g(e)  =  (v&w) 
(or  (v,w)  or  (w,v)  in  the  case  of  a  digraph),  then  e  is  said  to  be  in- 
cident to  v  and  w.   If  the  graph  is  directed,  then  an  edge  e  such  that 

g(e)  =  (v,w) 
is  said  to  be  incident  into  w  and  incident  out  of  v.   Alternately,  v  is 
referred  to  as  the  initial  vertex  of  e  and  w  is  referred  to  as  the  final 
vertex  of  e.   A  vertex  v  in  a  digraph  whose  incident  edges  are  incident 
out  of  v  is  called  a  source;  a  vertex  w  whose  incident  edges  are  incident 
into  w  is  called  a  sink. 

Geometrically,  we  represent  an  undirected  graph  by  associating  each 
vertex  with  a  geometrical  point  and  each  edge  with  a  curve  segment.   A 
digraph  is  represented  in  the  same  way,  except  that  the  curve  segment  is 
directed  according  to  the  order  of  the  vertices  corresponding  to  the  edge. 
Figure  1(a)  shows  a  geometrical  representation  of  an  undirected  graph, 
while  Figure  1(b)  shows  a  digraph. 


Figure  1(a)  an  undirected  graph 


(b)  a  digraph 


For    the   graph    in  Figure    1(a), 
V   - 


{WWV5J     • 
E  =  Ve^e^  e3,e4,e5,e6J      , 

g:e1    -      (v2&v3),    e2    -      (v3&v4) ,    e3    -      (v4&v5> ,    e4    -      (v^v^ , 

e5   -      (v^vp,    e6    -      (v1&v3). 

For   the   graph   in   1(b), 

V   -         {vj.v'.v'.v'.v'}     , 

E'    =  \-ei,e2,e3,e4,e5,eDj      ' 

g':    e|    -      (v{,v|),    e«    -      (vj,v^),    e'    -      (v^,v'),    e^   -      (vj.vj) 

e'    -      (v3',v2),    e^    -      (v',v|). 

In  the  sequel,  we  shall  consider  only  digraphs. 

In  a  digraph,  two  edges  are  said  to  be  consecutive  provided  the 

initial  vertex  of  one  is  the  final  vertex  of  the  other.   Thus,  two  edges 


e.  and  e9  such  that 


g(ex)  =  (a,b)  and  g(e2)  =  (b,c) 


are  consecutive  edges.   A  sequence  of  distinct  consecutive  edges  of  the 
form 


P  ■  (e1,e2,e3,...,en_1) 


such  that 


g(ex)  =  (v1,v2),  g(e2)  =  (v2,v3),  g(e3)  =  (v3»v4) . • • • >8(en_2^  = 

<Vn  9'v„  0      8(en  i>  =  (vn  i»v„) 
n-z  n-i      n-i      n-l  n 

is  called  a  path.   Note  that  no  edges  are  repeated  in  a  path,  however, 

vertices  may  be  repeated.   A  subsequence  P'  of  consecutive  edges  is 

called  a  subpath  of  the  path  P.   If  P'  starts  with  the  same  edge  e  as 

P,  we  say  P'  is  an  initial  subpath  of  P.   We  shall  say  a  vertex  v.  is 

attainable  from  a  vertex  v  provided  there  is  a  path  from  v  to  v.. 
o  r         o     i 
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Alternately,  we  say  that  v.  is  a  descendant  of  v  and  that  v   is  a  pre- 

i  o  o 

decessor  of  v  . 

We  define  a  function  d  on  the  set  ty   of  all  paths  in  a  graph  G  in- 
to the  set  of  natural  numbers  (including  zero)  by 

d(P)  =  n,  P  in  ff  , 
where  n  is  the  number  of  edges  in  the  path  P.   d(P)  is  referred  to  as 
the  length  of  the  path  P.   Note  that  an  edge  is  a  path  of  length  one. 
Let  G  =  G(V,E,g)  be  a  digraph.   Define  a  multi-valued  function  f 
on  V  into  V  according  to  the  following  rule: 

f(v)  =  i  v1:  v '  in  V  and  there  exists  e  in  E  such  that  g(e)  = 
(v,v')r   ,  for  all  v  in  V.   We  call  f  the  descendent  function  of  the 
vertex  v.   We  define  powers  of  f  as  follows: 


£J( 


v)  =  Jv':  v'  in  V,  v1  in  D   j  ,  v  in  V, 


where  D  .  is  the  set  of  all  descendants  of  v  by  a  path  of  length  j.   For 
example,  in  Figure  1(b),  we  have 

f(v{)  ■  [vj,  v'}  ,  £4(vJ)  =  [vj.v'.v-]  . 

A  special  type  of  path  in  a  digraph  is  one  which  starts  at  a  vertex 

v  and  ends  at  a  vertex  v  ,  where  v  =  v  .   Such  a  path  is  called  a  cycle 
o  n         o    n  — * 

If  no  vertex  (except  v  )  is  used  more  than  once  in  a  cycle,  then  the 

cycle  is  said  to  be  simple.  A  simple  cycle  of  length  one  is  called  a 

loop.   A  graph  which  contains  no  cycles  is  referred  to  as  acyclic. 

Suppose  G  »  G(V,E,g)  is  an  acyclic  graph  in  which  all  vertices  are 

attainable  from  a  single  vertex  v  in  V.   Further,  suppose  no  vertex  has 

more  than  one  edge  incident  into  it.   Then  G  is  called  a  rooted  ditree. 

v  is  called  the  root  of  G.   Figure  2  below  shows  a  rooted  ditree 
o  b 
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Figure  2.   A  rooted  ditree. 


Note  that  a  rooted  ditree  must  have  only  one  root.  A  vertex  in  a  rooted 

ditree  is  said  to  be  terminal  if  it  has  no  descendants.   Thus,  in  the 

graph  of  Figure  2,  a  is  a  root  and  a  ,a  ,a  ,a.  ,a. . ,a  are  terminal 

vertices.   A  branch  in  a  rooted  ditree  is  any  path  whose  initial  vertex 

is  the  root  and  whose  final  vertex  is  terminal.   Hence,  the  sequence  of 

edges 

B  =  (  (a0,a2),(a2,a5)(a5,a10)  ) 

(note  that  since  there  are  no  multiple  edges  between  any  pair  of  vertices, 

there  is  no  confusion  if  we  represent  an  edge  by  its  associated  pair  of 

vertices)  is  a  branch  of  the  ditree  in  Figure  2. 

A  graph  G  =  G(V,E,g)  whose  vertex  set  V  can  be  partitioned  into  two 

non-null  disjoint  sets  V,  and  V"  such  that  E  consists  only  of  edges  e. 
j  12  i 

such  that 


g(ei)  =  (vi»vj)  or  (v-tV^),  V±   in  V^  v  in  V2, 

is  called  a  bipatite  directed  graph  or  bi-digraph.   The  graph  in  Figure 
3  is  a  bi-digraph. 


Note  that 


V  -  [a1,a2,a3,a4,a5,a6,a7,a8,a9] 
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Figure  3.   A  bi-digraph 
We  have  the  following  result  concerning  bi-digraphs. 

THEOREM  I;   If  G  is  a  bi-digraph,  then  all  simple  cycles  in  G  are  of 
even  length. 

Proof:   Sup] ose  C  is  a  simple  cycle  in  G  of  the  form 

C  =  (e  ,e0,  ...,e   ,e  ), 
1  I  n-1  n 

where  d (C)  is  odd.   If  the  length  of  C  is  odd,  then  the  initial 

vertex  of  C  must  be  different  from  the  final  vertex  by  the  definition 

of  bi-partite  graph.   But  this  is  impossible  since  the  initial  and  final 

vertex  of  a  cycle  must  coincide.   Therefore,  d (C)  is  even. 
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III.   AN  ALGORITHM  FOR  FINDING  THE  SIMPLE  CYCLES  IN  A  DIGRAPH 

An  important  problem,  especially  in  the  application  of  graph  techni- 
ques, is  that  of  locating  simple  cycles  in  a  digraph.   In  dealing  with 
large  graphs  (large  in  the  sense  that  V  and  E  contain  many  elements), 
it  is  useful  to  have  an  efficient  algorithm  which  conserves  memory  space 
of  a  computer.   The  algorithm  we  shall  develope  in  this  chapter  is  de- 
signed to  meet  these  requirements.   In  the  final  chapter  of  the  paper,  we 
will  show  how  our  algorithm  compares  with  an  algorithm  given  by  Denielson 
[3]. 

Before  starting,  it  will  be  necessary  to  introduce  some  new  terminol- 
ogy.  In  the  sequel,  we  shall  denote  the  vertices  of  all  graphs  under 
consideration  by  intef^r-indexed  lower  case  letters  such  that  no  two 

vertices  have  the  same  index.   If  j  is  an  integer  such  that  v.  occurs  in 

J 

a  graph  under  consideration,  we  will  say  v.  is  the  vertex  corresponding 

to  the  indexing  integer  j. 

•. 

DEFINITION  1. 

Let  G  =  G(V,E,g)  be  a  digraph.   Let  v.  be  an  element  of  V.   Then, 
the  indexing  integer  j  corresponding  to  v .  is  said  to  be  an  adjacent 
element  of  v  provided  v.  belongs  to  f(v  ),  where  f  is  the  descendent 
function  of  v  . 

For  example,  in  the  graph  in  Figure  3  above,  the  adjacent  elements 
of  vs  are  the  integers  6  and  9. 

DEFINITION  2. 

If  G  =  G(V,E,g)  is  a  digraph,  and  if  v.  is  a  vertex  in  V,  then  the 
adjacency  set  of  v  is  the  set  of  all  adjacent  elements  of  v.  (e.g.,  the 
adjacency  set  of  v  is  the  set  of  index-integers  corresponding  to  the 
elements  of  f(v  )). 


For  example,  in  the  graph  in  Figure  3,  the  adjacency  set  of  a  is 
the  set  (6,9j  ■  If  no  vertices  are  adjacent  to  a  vertex  (i.e.,  if  a  ver- 
tex is  a  sink)  we  denote  its  adjacency  set  by  0  (the  null  set). 

In  the  discussion  which  follows,  we  assume  that  the  graph  under  con- 
sideration, G  =  G(V,E,g)  is  a  digraph  on  n  vertices  which  are  represented 
by  the  lower  case  letter  v  indexed  1  through  n.   We  further  assume  that 
in  order  to  communicate  to  the  computer  the  structure  of  G,  we  have  a 
list  which  associates  each  vertex  in  V  with  its  adjacency  set. 

DEFINITION  3. 

An  attainability  chain  of  a  vertex  v  in  a  digraph  G  is  a  directed 
graph  C. ,  which  corresponds  to  a  simple  path  P(v . ,v.)  in  G.   The  source 
vertex  of  C.  is  i  (corresponding  to  v  ),  and  the  sink  vertex  is  either  a 
sink  vertex  in  G,  or  a  copy  k'  of  the  first  repeating  vertex  v  in  the 
path. 

For  example,  consider  the  graph  in  Figure  3.   The  vertex  a  has  the 
following  attainability  chains: 

im  1 

3<L       3 

4*      4 

56       2 

6  0       1  "< 

Figure  4.   Attainability  chains  for  the  vertex  a  of  the  graph  in  Figure  3 
Note  that  if  a  simple  cycle  occurs,  as  in  the  cycle  composed  of  the  ver- 
tices a.,  a„,  a,,  a  ,  a  ,  then  the  attainability  chain  of  any  of  these 
vertices  terminates  as  soon  as  the  first  repetition  of  a  vertex  occurs. 
Also,  note  that  each  attainability  chain  is  a  rooted  ditree  with  one 
branch. 
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DEFINITION  4. 

The  attainability  tree  for  a  vertex  v.  is  a  rooted  ditree  wi'h  a 
vertex  corresponding  to  i  in  the  root  position  and  with  branches  con- 
sisting of  all  the  attainability  chains  of  v  .  We  will  super-impose  the 
initial  subpaths  of  any  attainability  chains  with  identical  initial  sub- 
paths. 

The  graph  in  Figure  5  shows  the  attainability  tree  for  the  vertex  a. 
of  Figure  3. 
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Figure  5.  An  Attainability  Tree 

We  now  give  the  formal  statement  of  the  algorithm,  keeping  in  mind 
that  the  graph  under  consideration  is  a  digraph  G  =  G(V,E,g)  on  n  ver- 
tices, and  that  we  have  a  list  of  vertices  and  their  corresponding 
attainability  sets. 

STEP  1,    Inspect  the  list  and  eliminate  all  vertices  whose  index 
integers  do  not  appear  as  adjacent  elements.   Obviously,  these 
vertices  could  not  belong  to  simple  cycles  because  they  have  no 
edges  incident  into  them  (i.e.,  they  are  sources).   Also,  inspect 
the  list  and  eliminate  as  roots  all  vertices  whose  adjacency  set  is 
0.   Again,  it  is  clear  that  such  vertices  caanot  belong  to  any 
simple  cycles  because  there  are  no  edges  incident  out  of  them 
(i.e.,  they  are  sinks).   We  have  thus  eliminated  all  sourcas  and 
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sinks  as  possible  roots.   We  may  now  repeat  the  procedure  on  the 
remaining  vertices  until  no  more  eliminations  are  possible. 

STEP  2.   Pick  any  vertex  v   from  the  list  which  has  not  been  eli- 
minated and  construct  its  attainability  tree,  consulting  the 
origional  list  for  incidence  relations.   Recall  that  we  super- 
impose identical  initial  subpaths. 

STEP  3.   If  an  element  is  repeated  along  any  attainability  chain, 
record  that  element  together  with  the  elements  between  the  two 
occurrence  (in  the  order  they  appear).   The  vertices  corresponding 
to  these  elements  will  compose  a  simple  cycle  in  G.   Note  that 
repeated  vertices  must  occur  along  the  same  branch  of  the  at- 
tainability tree  for  a  cycle  to  be  located. 

STEP  4.   Eliminate  as  possible  roots  all  elements  which  appeared  in 
the  attainability  tree  of  v  ,  for  if  these  vertices  belonged  to 
simple  cycles,  this  would  have  been  discovered  in  the  attainability 
tree  of  v  . 

STEP  5.   If  all  vertices  have  been  eliminated  as  possible  roots, 

the  algorithm  terminates.   If  some  vertices  remain,  choose  any  such 

vertex  v,  not  eliminated  as  a  possible  root,  construct  its  attain- 
k 

ability  tree  and  return  to  step  3. 

We  now  ask  the  following  question:   Does  the  algorithm  locate  all 
the  simple  cycles  in  a  digraph  (and  only  the  simple  cycles)?  We  answer 
this  question  with  the  following  theorem: 
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THEOREM  II.  The  algorithm  outlined  in  Steps  1  through  5  above  locates  all 
of  the  simple  cycles  of  a  digraph,  and  it  locates  only  simple  cycles. 

Proof :  First,  note  that  no  cycles  will  be  located  which  are  not 
simple,  since  a  branch  is  terminated  at  the  first  repetition  of  a  vertex. 

Now,  suppose  there  is  a  simple  cycle  C  in  the  digraph  G  which  was 
not  located  by  the  algorithm.   Then  C  can  have  no  edges  in  common  with 
any  attainability  tree  considered  in  the  sequence,  for  if  it  did,  then 
by  the  fact  that  a  branch  is  terminated  only  if  a  source,  a  sink,  or  a 
repetition  occurs,  C  would  have  been  located.   Conversely,  only  those 
edges  of  G  which  involve  sources  or  sinks  in  G  are  left  out  of  the  set 
of  attainability  trees.   Thus,  we  have  established  that  if  C  is  not 
located  by  the  algorithm,  then  C  does  not  belong  to  G,  which  is  a  contra- 
diction.  Therefore,  C  must  have  been  located  by  the  algorithm. 

It  should  be  noted  at  this  point  that  some  cycles  may  be  located 
more  than  once.   However,  this  will  occur  only  in  special  cases. 

As  an  example  of  the  application  of  this  algorithm,  consider  the 
digraph  shown  in  Figure  6  below. 


Figure  6.   A  Digraph 
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We  have    the    following   list  of  vertices  with   their  adjacency   sets: 


v  f105 
v  t11} 
v  fu) 

v?:   [l3,    «} 
V   {w} 

w 

M 

to 

0 


10 


11 


12 


13 


Applying  Step  1,  we  eliminate  v.,  as  a  possible  root,  since  its  adjacency 
set  is  0.   We  also  eliminate  v.,  v  ,  v  ,  and  v.,  since  1,  2,  3,  and  4 
do  not  occur  in  any  of  the  adjacency  sets.   Let  us  now  construct  the 


attainability  tree  of  v 


SIN  K 

;tainability    tree    for  vin  of  Figure   6 


10 
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We  have  located  one  simple  cycle,  namely 

Note  that  5,  6,  7,  8,  9,  10,  11,  12,  13  all  appeared  in  the  attainability 
tree  of  v   ,  and  that  1,  2,  3,  4,  and  14  have  already  been  eliminated 
as  roots.   Hence,  by  Step  4,  we  need  construct  no  more  attainability 
trees,  and  we  can  conclude  that  C  is  the  only  simple  cycle  in  the  graph. 
An  inspection  of  Figure  6  also  leads  to  this  conclusion. 

Note  that  if  the  graph  under  consideration  were  a  bipartite  digraph 
(or  bi-digraph) ,  we  would  have  had  two  lists  instead  of  one.   In  this 
case,  we  would  add  the  following  step  the  simplifications  above: 

STEP  0:   Eliminate  as  possible  roots  all  vertices  which  occur  in  the 
list  with  the  most  vertices.   This  will  be  sufficient  to  locate  all 
simple  cycles  by  Theorem  I  which  states  that  all  cycles  in  a  bi- 
digrapl  have  even  length,  and  hence  must  contain  at  least  one  vertex 
from  each  list.   This  simplification  will  be  used  in  the  next 
chapter. 

We  observe  at  this  point,  that  though  the  primary  purpose  of  our 
algorithm  as  stated  is  to  locate  simple  cycles,  it  may  also  be  used  to 
locate  any  simple  path  in  a  digraph.   For  example,  if  we  wish  to  know 
if  there  is  a  simple  path  between  a  vert-x  v0  and  a  vertex  v.  in  a  di- 
graph, we  need  only  construct  the  attainability  tree  of  vQ  and  note 

whether  it  contains  v.. 

l 
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IV.   COALESCED  GRAPHS 

The  purpose  of  the  foregoing  algorithm  is  to  make  possible  the  loca- 
tion of  all  simple  cycles  in  a  digraph  with  a  large  number  of  vertices. 
In  practice,  in  the  analysis  of  very  large  graphs,  such  as  the  graph  of 
the  circuitry  of  a  digital  computer,  a  procedure  is  used  which  partitions 
the  graph  into  subgraphs  which  are,  in  turn,  partitioned  into  smaller  sub- 
graphs, until  manageable  subgraphs  are  obtained.   Once  these  subgraphs 
have  been  analyzed  for  cyclic  structure,  they  are  generally  coalesced 
into  a  single  vertex,  and  the  resulting  graph  is  then  analyzed  for  its 
cyclic  structure.   In  this  chapter,  we  will  show  how  the  technique  of 
Chapter  II  may  be  applied  to  such  "coalesced  graphs." 

One  of  the  principal  problems  in  locating  the  simple  cycles  in  a 
coalesced  graph  is  that  a  vertex  which  is  a  coalesced  subgraph  may  be 
attainable  from  a  vertex  v,  but  a  vertex  v*  which  is  attainable  from 
the  coalesced  subgraph  may  not  be  attainable  from  v  in  the  original 
graph.   For  example,  suppose  G  is  the  graph  shown  in  Figure  7. 


Figure  7.  A  Digraph 
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Suppose  that,  by  some  procedure,  we  determine  that  the  subgraph  of  G 
composed  of  the  vertices  a,  b,  c,  d,  e,  and  f,  together  with  the  edges 
connecting  them,  can  be  coalesced  into  a  single  vertex  v  .   Then,  the 
graph  shown  in  Figure  8  is  the  coalesced  version  of  the  graph  in  Figure 
7. 


-® Ht 


V* 


Figure  8.   The  Coalesced  Version  of  the  Graph  in  Figure  7 
Note  that  this  graph  is  bipartite,  since  one  node  actually  represents  a 
subgraph,  while  the  others  represent  single  vertices  in  the  original 
graph.   Also,  note  that  in  the  graph  in  Figure  8,  v*  is  a  descendant  of 
v  but,  in  the  original  graph  (Figure  7),  it  is  not.   Thus,  we  see  that 
our  algorithm  cannot  be  applied  to  coalesced  graphs  without  some  modi- 
fication. 

Let  G  =  G  (V  ,E  ,g  )  be  a  digraph  without  coalesced  vertices, 
o    o  o  o,&o        ° 

Let  S  be  a  subgraph  of  G  .   Form  the  graph  G  ■  G  (V  ,E  ,g  )  in  which 

o  c    c  c  c  c 

S  has  been  coalesced  into  a  single  vertex. s.   Let  f  be  the  descendent 

o 

function  on  the  vertices  of  G  ,  and  let  f  be  the  descendent  function  on 

o  c 

the  vertices  of  G  .   Define  subsets  I(S)  and  J(S)  on  the  vertex  set  V 
c  c 


of  G  as  follows: 

c 


KS)  =  £Vj,  in  Vc:  s  is  in  f<.(vi)J 

J(S)  =   Vv,  in  V  :  v.  is  in  f  (s)| 
»-  j     c   j       c   J 


20 


Now,  define  a  mapping  p  on  I(S)  x  J(S)  into   0,1  such  that 


0  if  v.  is  not  in  fk(vj  for  all  k, 

1  if  there  exists  k  such  that  v.  is  in  f  (v.) 

J        °   i' 


Compute  p(v  »v.)  for  all  pairs  (v., v.)  with  v  in  I(S)  and  v .  in  J 
■*■   J  •*■   J        ■*•  J 


(S) 


as  follows:   Let  s.  be  in  S  such  that  s.  is  also  in  f  (v.).   Let  s.  be 

l  l  o  i         j 

such  that  v.  is  in  f  (s.).   Note  that  s.  and  s.  must  exist  since  v.  is 
J        o  j'  i      j  l 

in  I(S)  and  v.  is  in  J(S).   Let  T   be  the  attainability  tree  of  s  , 

i 

and  suppose  V   is  the  vertex  set  of  T   .   Then,  if  s.  is  in  V   , 
r     s.  s.  j       s. 

l  i  J        l 


otherwise, 


P(v..,v  )  -  1, 


P(v.v .)  =  0. 

i  y 


Note  that  T   may  not  exist,  in  which  case  p(v_,v.)  is  obviously  zero. 
s.  i  j 

l  J 

When  p(v.,v.)  has  been  computed  for  all  ordered  pairs  (v., v.)  with 
l  j  i  j 

v  in  I(S)  and  v.  in  J(S) ,  we  apply  the  algorithm  of  Chapter  II  to  the 

bi-digraph  G  .   Suppose  we  do  this  and  obtain  simple  cycles  C  ,  ...  ,  C  . 

Note  that  by  Theorem  I,  the  length  of  C.  (i  =  1,  ...  ,  n)  will  be  even. 

Hence,  if  nu  is  the  number  of  vertices  in  C. ,  then  m./2  is  the  number  of 
i  11 

coalesced  subgraphs  in  C..   Denote  these  coalesced  subgraphs  by  k.,,kJ_, 

l  °  r  ll   i2 

"    *    kimi/2* 
Let 

I(k.  .)  =  )v   in  V   :  k.  .  is  in  f  (v  )f  ,  i  =  1,  ...  ,  n;  j  =  1,  . ..,m./2 
ij    I  m     c    ij        c  n  J  J  i 

J(k.  .)  =  tv.  in  V   :  v.  is  in  f  (k. .)j  ,  i  =  1,  . . .  ,  n;  j  =  1,  ...,m./2 
ij    i   1     c    1        c   ij  J  '      ■     '    J  '  '  1 

For  each  k. .  (i  =  1,  ...,  n;  j  =  1,  ...,  m  /2)  form  the  ordered  pair 

(t .  - ,t. .-),  where  t .  .  is  in  I(k. .)  and  t . ,,  is  in  J(k. .)  and  both  t .  , 
J-l   J+l  J-l         ij       j+1  ij  J-l 

and  t    are  in  C  .   Now,  the  ordered  pair  (t   .,  t   )  must  be  in  the 
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domain  of  p,  since  (t.^.t.^)  is  in  1(1^  )  x  J(k_)  .   Note  the  value  of 
p(t.   ,t   ).   If  it  is  zero  for  any  j  (j  -  1,  ...  ,  mj2) ,    then  (^  is 
not  actually  a  cycle  in  G  .   If  it  is  one  for  all  j,  then  C,.  is  indeed  a 


cycle  in  G  .   For  example,  suppose  G  is  the  graph  in  Figure  9  below. 

J  r\  O 


\Z*.n 


Vx9    \/st 


V4t 


VCT     Vc3 


Vt( 


'7r 


Figure  9.   A  Digraph 
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This  is  a  digraph  of  an  adder  from  a  digital  computer.   Suppose  we  have 
a  procedure  which  partitions  Gq  into  the  four  subgraphs  shown  in  Figure 


tfu 


«n  v*r 


v^ 


V*1       SUG&RftPK  (?*  ^Jf 


v£t   sua  £*aph   c 


Figure    10. 


Four  Subgraphs   of   the  Graph    in  Figure   9. 
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We  now  apply  our  algorithm  to  the  four  subgraphs.   We  have  the  following 
lists  of  vertices  together  with  their  adjacency  sets: 


SubR 

raph  A 

Subgraph   B 

Subg 

raph  C 

SubR 

raph   D 

V2   : 

M 

v!    ••   i27} 

V15: 

£53,54,553 

V8    : 

125^ 

V3    : 

M 

v8    1   {24} 

v16: 

{53,543 

V9    : 

{25} 

V4    : 

M 

v9    1   [24} 

v17: 

{66,68/ 

Vll: 

{79} 

V5    : 

M 

V!0=    W 

v18: 

{54} 

V12: 

{79} 

V6    : 

hi 

v27:  fc8? 

V19: 

{33} 

V13: 

{80} 

V7    : 

[23} 

v24:  {36.37} 

v      : 

20 

{56} 

V14: 

{80} 

V21: 

{»} 

v28:  M 

V53: 

{59,653 

V25: 

{38,39} 

V22; 

[29] 

v36:  [32,33] 

V54: 

{60} 

V79: 

{26) 

v23: 

{29} 

v35:  H 

V55: 

M 

v      : 

80 

{26} 

V29: 

{30} 

v40:   C33} 

V56: 

M 

V35: 

{76} 

V30: 

{35} 

v37:   W 

v63: 

{.66,67} 

v      t 

38 

[76,77) 

V35: 

0 

v32=  [«,«} 

V59: 

{68} 

V39: 

£78} 

V63: 

H 

v33:  {41,43} 

V60: 

}66,67,68,< 

>9}  V26: 

f3l] 

V71: 

{A 

v34:    {41,44} 

V65: 

^67,69*} 

V31: 

<Uo} 

V72: 

M 

v42=    H 

V64: 

{68,69} 

V40: 

{773 

V70: 

H 

v43:   W 

V66: 

{75} 

vio: 

{.78} 

v?3: 

H 

v44:    {so} 

V67: 

{"} 

V76: 

{45,46} 

V74: 

H 

v50:  [51> 

v68; 

M 

V77: 

{45,473 

V75: 

0 

\v  hi 

V69: 

M 

V78: 

[45,48} 

V61: 

{«} 

v45=    &7} 

V75: 

0 

V46: 

{52l 

V62: 

t«l 

v51:  {34,58} 

V71: 

0 

V47: 

£52} 

V81: 

M 

v5g:   {62} 

V72: 

0 

V48: 

t«i 

V64: 

0 

v5?:    {61,633 
v61,v62,v63:  ' 

V49: 
V52: 

{78} 
{49} 
0 
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We  have  the  following  attainability  treer,  for  the  subgraphs  A  through  D: 

•M        mXZ     *Xi      -C3 


EHO 


1  I  fl         ?*"'  ?"""         tU 

_w  r °  -i*  i»<  in-  L 


3^  4o    4-r 


34^J3  «      0  ST 


OfcLe:^,  *»,**/*) 


Atti^ih aQicvtv    t«.s&S     Foa   \5        


end      e*o 


KND 


•53 


Ac  A 


21 


It 

(UIO      RMO 


ore-*-,       tin  • 


Figure   11.      Attainability  Trees    for   the   Subgraphs   of  Figure    10. 
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We  have  located  two  simple  cycles: 

Cl  =  (V50'V51'V34'V44'*50) 

C2  =  (V78'V48'V52»V49'V78) 
We  now  coalesce  the  four  subgraphs  into  vertices  and  examine  the  re- 
sulting bi-digraph  for  cyclic  structure  in  order  to  discover  if  there  are 
any  larger  simple  cycles  in  G  .   Figure  12  shows  the  coalesced  graph  G  . 


Figure  12.  The  Coalesced  Graph  of  the  Graph  in  Figure  9. 
Since  this  is  a  bipartite  graph,  we  have  two  lists  of  vertices: 


A:  £35,64,751 

B:  £41,62,61,63i 

C:  {75,71,72} 

D:  £40,45} 


11 


M 

M 

La) 

W 

to 
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12 

r13 

'l4: 

r15 

r16 

r17 

r18 

'19 

r20 


'41 
'62 : 
r61: 

'35 

'75 


Cc} 

M 

fc} 
W 

w 

M 

0 


71 
'72 
r64; 
r63 


'40 


'10 
'45 ; 


M 

Ja,c} 
Vb} 

Cb,d3 

\b,dJ 

{b5 


We  now  invoke  step  0  and  consider  as  possible  roots  only  those  vertices 


which  appear  in  list  V  .   We  construct  the  attainability  tree  for  vertex 
A. 


Figure  13.   Attainability  Tree  for  Vertex  A   of  G 


CVCcR   OCCK 
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We  see  that  vertices  B,  C,  and  D,  appear  in  the  attainability  tree  of  A, 
so  it  is  unnecessary  to  construct  any  more  attainability  trees  (by  Step 

4). 

We  have  located  the  following  simple  cycles  in  G  : 

c 

,vf2,A) 

r61'A) 

'V63'A> 
(v63,C,v71,A) 


Cl  =<A>V35>B"62: 
C2  =(A,v35,B 

C3  =(A,v35,B 


=  (A,v.t.,B 


C5   =(A,v35,B 
C6  =(A»V35>D 

C7  =(A>V35'D 
C8  =(A,v35,D 

C9  =(A,v35,D 

C10  =(A'V35'D 

Cll  =(A>V35>D 
C12  =(A,v35,D 

C13  =(A,v35,D 

C14  =(A,v35,D 

C15  =(A'V35»D 

C16  =(A'V64>C 
C1?  =(A,v64,C 


v63,C,v?2,A) 

V40'B»V62'A> 

V40'B'V61'A) 

V40>B'V63'A> 

V40'B'V63'C'V71'A) 

V40'B'V63»C»V72»A> 

V45,BjV62'A) 

V45'B'V61'A) 
v45,B,v63,A) 

V45'B'V63'C»V71'A) 

V45»B'V63'C'V72'A> 
v71,A) 

v?2,A) 


We  must  now  compute  the  values  of  p  to  see  which  of  these  simple  cycles 

in  G  is  also  a  simple  cycle  in  G  .   We  shall  do  this  only  for  C, .   The 
c  o  1 

others  would  be  done  exactly  the  same. 

We  first  note  that  C   involves  the  coalesced  vertices  A  and  B,  with 
A  in  f  (v62  ),  v    in  f  (A),  B  in  f  (v  oc.),  and  vto  in  f„(B).   Hence, 


35' 


62 
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we  must  compute  P(V62»V35)  andp^v35»V62^ '   T°  comPute  P^V62,V35^'  We 
examine  subgraph  A,  and  note  that  62  does  not  appear  in  any  of  the  at- 
tainability trees  of  the  vertices  of  A.   Hence,  we  must  construct  the 

attainability  tree  of  v,„  to  see  if  v._  is  a  descendant  of  v.„: 

62  35  62 


Figure  14.  The  Attainability  Tree  of  v,_ 

We  see  that  35  does  not  appear  in  the  attainability  tree  of  v,  and  thus, 

p(vfi_,v  )  =  0.   Therefore,  C.  is  not  a  cycle  in  G  .   We  will  not  repeat 

the  procedure  for  C_  through  C.-,  however,  if  we  did,  we  would  find  that 

none  of  these  are  cycles  in  G  . 

o 

In  actual  applications  of  this  technique,  we  would  first  determine 

the  cyclic  structure  of  each  subgraph  which  is  coalesced  into  a  node. 

In  addition  to  retaining  the  resulting  information  about  cycles,  we 

would  also  note  p(v..v.)  for  each  pair  of  nodes  external  to  the  coalesced 

i  j 

subgraphs.   We  would  then  use  this  information  in  constructing  the  attain- 
ability trees  for  the  coalesced  graph.   Had  we  done  this  above,  the  at- 
tainability tree  in  Figure  13  would  have  been  less  complicated.   We  did 
not  folio. /  this  procedure  above  in  order  to  better  illustrate  the  use  of 
the  function  p. 
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V.   A  COMPARISON  WITH  ANOTHER  ALGOTHITHM 

We  shall  now  show  how  the  algorithm  developed  in  Chapter  II  for 
locating  the  simple  cycles  of  a  digraph  is  much  superior,  in  terms  of 
operations  required,  to  another  algorithm  used  for  the  same  purpose. 

Danielson  [3]  developed  an  algorithm  for  locating  all  paths  (in- 
cluding cycles)  in  an  undirected  graph  using  matrix  techniques.   His 
method,  while  designed  for  undirected  graphs,  can  easily  be  modified 
for  application  to  digraphs.   We  now  state  the  necessary  definitions  to 
explain  his  algorithm. 

DEFINITION  5. 


If  G  -  G(V,E,g)  is  an  undirected  graph  on  n  vertices  (v  ,  ...  ,  v  ) 

1        n 

and  m  edges  (e  ,  ... ,e  ) ,  then  the  adjacency  matrix  A  of  G  is  given  by 


A  =  (a. .) ,  where 
o     ij 


au  = 


1  if  an  edge  connects  v.  to  v , , 
0  otherwise. 


DEFINITION  6. 


If  G  =  G(V,E,g)  is  an  undirected  graph  on  n  vertices  (v.,  ...  ,  v  ) 

1        n 

and  m  edges  (e  ,  ...  ,  e  ) ,  then  the  variable  adjacency  matrix  A  of  G 


is  given  by  A  =  (a. .),  where 


ij 


(ekif 
^0   o 


edge  e,  joins  v.  to  v., 
k        i     j' 

therwis  . 


DEFINITION  7. 


If  G  =  G(V,E,g)  is  an  undirected  graph  on  n  vertices  (v  ,  .. 
then  the  modified  variable  adjacency  matrix  B  of  G  is  given  by  B 
where 


•V- 

<»  >, 
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ij 


[v.   if  an  edge  joins  v,  to  v  . , 
I  0   otherwise 


To  illustrate  these  concepts,  consider  the  graph  of  Figure  15. 


We  have 


Figure  15.  An  Undirected  Graph 


A  = 
o 


0  10  0  1 
10  0  11 
0  0  0  11 
0  110  0 
1110  0 


,  A  = 


0     e,    0 

0 

e. 

0 

v„ 

0 

0 

vc 

1 

2 

2 

5 

e,    0     0 

e^ 

e„ 

V, 

0 

0 

v. 

v_ 

1 

6 

3 

1 

4 

5 

0     0     0 

e5 

e4 

,    B   = 

0 

V2 

V3 

0 

0 

0     e     en 

0 

0 

0 

v„ 

v^ 

0 

0 

6     5 

2 

3 

e„   e„   e. 

0 

0 

vn 

v. 

v. 

0 

0 

-2     3     4 

w. 

_  1 

2 

3 

The  algorithm  states  that  all  edge  progressions  of  length  k  between 
v  and  v.  will  appear  in  the  i-j  position  of  P  ,  where  P  is  defined  re- 
crusively  as  follows: 


P  =  BA  .  P  =  BP,   . 
2     ok     k-1 
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Hence,  for  the  graph  in  Figure  15, 


VS^5V2 

$  1 

Vila's 

V2V3 

V1 

TuWz 

V5Y1 

VS 

TiWW3 

T,T2V2Y3 

V 

TuW2 

V1 

0 

V2V3 

Y2V3 

T2Ws 

V^3V5 

VuVs 

0 

V2V1 

V5V5 

Vs^Ws 

Vg^aVz's 

V« 

T,T2  ^Z'l 

We  first  note  that,  as  stated  before,  the  algorithm  only  locates 
the  paths  and  cycles  in  an  undirected  graph.   However,  it  is  possible  to 
modify  the  procedure  slightly  to  obtain  the  cycles  of  a  directed  graph. 
How  this  is  done  is  unimportant;  what  is  important  is  that  to  locate  all 
paths  (including  cycles)  of  lengths,  k,  k-1,  ...  ,  2  in  a  digraph,  we 
are  required  to  multiply  n  x  n  matrices  k  times.   In  particular,  if  we 
were  seeking  all  simple  cycles  in  a  complicated  graph  on  n  vertices, 
we  would  have  no  choice  but  to  perform  matrix  multiplication  n  times. 
Now,  in  using  this  algorithm,  we  are  not  performing  matrix  multiplication 
in  the  usual  manner,  since  the  elements  of  the  P,  's  and  of  B  are  symbols 
(i.e.,  subscripted  lower  case  letters)  rather  than  numbers.   Hence,  what 
we  are  really  doing  is  concatenating  symbols  in  order  to  obtain  the  ele- 
ments of  the  P  's.   Note  that  the  symbol  +  is  included  in  this  concaten- 
ation each  time  strings  of  symbols  are  added  together.   For  example,  the 
element  of  P_  above  which  appears  in  the  second  row  and  fourth  column 
(VjV_-+v  v.+v  v_)  is  the  concatenation  of  eight  symbols.   We  interpret 
this  string  of  symbols  as  follows:   there  are  three  paths  of  length  three 
from  vertex  v_  to  vertex  v, ,  namely  via  vertices  v  and  v_  or  via  vertices 
v,.  and  v_  or  via  v  and  v,,.  We  observe  that  the  basic  operation  of 
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Danielson's  algorithm  is  the  writing  down  of  a  single  symbol.   That  is, 
in  writing  the  element  mentioned  above,  eight  basic  operations  were  in- 
volved.  This  process  takes  up  more  time  and  memory  space  than  that  re- 
quired to  perform  a  comparison  of  two  integers  which  is  the  basic  operation 
of  our  algorithm. 

As  was  mentioned  earlier,  we  are  required  to  perform  n  nxn  matrix 

multiplications  in  order  to  analyze  an  arbitrary  graph  for  cyclic  struc- 

3 
ture  using  Danielson's  algorithm.   Hence,  we  must  generate  n  elements, 

each  of  which  may  be  a  long  string  of  symbols.   Though  some  of  these  ele- 
ments may  be  zero,  at  worst,  as  in  the  case  where  each  pair  of  vertices 
is  joined  by  an  edge  and  a  loop  occurs  at  each  vertex,  all  elements  would 

be  non-zero.   Now,  in  the  event  we  were  to  analyze  such  a  graph,  each 

2  k-1 

element  of  P  would  contain  n  terms  of  (k-l)n   '  symbols  for  the  nodes 

k-1  2 

and  n   '  -  1  plus  (+)  signs.  Altogether,  each  of  the  n  terms  of  P 

K 


would  require 


(k-l)nk'1+  hk"L-l  =  knk-1  -  1 


symbols.   Hence,  if  the  basic  operation  is  that  of  writing  down  a  single 
symbol,  then  the  number  of  operations  required  to  locate  all  paths  (and 
thus,  cycles)  in  the  worst  case  is  given  by 


D  -  J^n1"1  -  l)  . 


k=2 

To  obtain  an  upper  bound  for  the  number  of  basic  operations  (i.e., 
comparison  of  integers)  using  our  algorithm,  we  consider  a  special  type 
of  graph,  namely  the  complete  symmetric  graph  on  n  vertices.  A  digraph 
G  =  G(V,E,g)  is  said  to  be  complete  if  for  every  pair  of  vertices  v.  and 
v  in  V,  there  exists  an  edge  e  in  E  such  that  either  g(e)  =  (v  ,v.),  or 
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g(e)  =  (v.,v  ).   A  digraph  G  =  G(V,E,g)  is  said  to  be  symmetric  if 
for  every  edge  e  such  that  g(e)  =  (v., v.)  for  v  and  v.  in  V,  there 
exists  an  edge  e'  such  that  g(e')  =  (v  ,v.).   Thus,  a  complete  sym- 
metric graph  on  4  vertices  is  shown  in  Figure  16  below: 


Figure  16.  A  Complete  Symmetric  Graph  on  4  Vertices 
We  will  compute  the  number  of  comparisons  N  necessary  to  analyze  this 
graph  for  cycles.   Clearly  this  will  be  an  upper  bound  on  the  number  of 
comparisons  required  to  analyze  any  digraph  without  parallel  edges,  since 
any  other  digraph  with  these  properties  will  be  a  subgraph  of  the  complete 
symmetric  graph  on  the  same  number  of  vertices. 

Note  that  the  adjacency  set  of  each  vertex  in  a  complete  symmetric 
graph  will  contain  n  elements.   Hence,  at  the  first  stage  of  the  attain- 
ability tree  (i.e.,  at  the  stage  involving  all  vertices  attainable  by  a 
path  of  length  one  from  the  root  vertex),  we  have  n  elements  which  must 
be  compared  with  one  element,  namely  the  root  vertex.  At  this  stage, 
we  will  locate  one  simple  cycle:   the  loop  which  occurs  at  the  root  vertex, 
Thus,  we  continue  on  to  the  second  stage  with  n(n-l)  vertices.   But,  at 
this  second  stage,  we  will  locate  all  of  the  other  n-1  loops  and  all  of 
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the  simple  cycles  of  length  two  which  involve  the  root  vertex.   There 
are  n-1  of  these.   Therefore,  we  continue  on  to  the  third  stage  with 
(n(n-l)  -  (n-1)  -  (n-1))  n  =  (n-1)  (n-2)n  vertices.   Now,  at  this  stage, 
we  locate  (n-1)  (n-2)  loops  (obviously  some  are  duplications),  (n-1)  (n-2) 
cycles  of  length  two,  and  (n-1)  (n-2)  cycles  of  length  three.   Thus,  we 
continue  on  to  the  fourth  stage  with 

(n(n-l)(n-2)  -  (n-1) (n-2)  -  (n-1) (n-2)  -  (n-l)(n-2))n 
=  (n(n-l)(n-2)  -  3(n-l) (n-2))n  =  (n-1) (n-2) (n-3)n 
vertices.  At  the  kth  stage,  we  will  have  (n-1) (n-2)  ...  (n-k+l)n  vertices. 
Now,  each  of  these  must  be  compared  with  k  vertices.   Hence,  at  the  kth 
stage,  there  are  kn(n-l)(n-2)  ...  (n-k+1)  comparisons  to  make.   We  demon- 
strate this  more  formally  in  the  following  theorem. 

THEOREM  III.   For  a  complete  symmetric  digraph  (with  loops  at  each  vertex) 
on  n  vertices,  there  will  be  (n-1) (n-2)  ...  (n-k+l)n  vertices  at  the  kth 
stage  of  the  attainability  tree. 

Proof:   The  proof  is  by  induction  on  k.   Suppose  k  is  1.   Then  the 
number  of  vertices  is  n,  which  is  clearly  the  case,  since  in  the  complete 
symmetric  graph  on  n  vertices,  each  element  has  n  terms  in  its  adjacency 
set. 

Now,  suppose  the  conclusion  holds  for  k  =  p.   That  is,  suppose  that 
at  the  pth  stage,  there  are  (n-1) (n-2)  . ..(n-p+l)n  vertices.  At  this  stage, 
we  will  locate  (n-1) (n-2)  ...  (n-p+1)  loops,  a  like  number  of  cycles  of 
lengths  2,  3,  ...  ,  p.   Hence,  at  the  p+1  st  stage,  we  will  have 

((n(n-l)(n-2)  ...  (n-p+1)  -  p(n-l)(n-2)  ...  (n-p+l))n 
=  n(n-l)(n-2)  ...  (n-p+1) (n-p) 
vertices.   Thus,  since  p  was  arbitrary,  the  conclusion  holds  in  general. 
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We  see,  then,  that  an  upper  bound  on  the  number  of  comparisons  re- 
quired to  analyze  an  n-vertex  digraph  is  given  by 


=  2il   kn(n-l).7(n-k).' 


C 

k=l 

k-1 
Now,  since  n   "  is  clearly  greater  than  (n-1)  .'/(n-k) .'  for  k  greater  than 

one,  C  will  be  less  than  D  (for  n  "5:  3) .   In  fact,  as  n  becomes  large, 
the  difference  between  D  and  C  becomes  very  great  indeed.   Also,  since 
the  complete  symmetric  graph  occurs  rarely  in  practice  the  number  of 
comparisons  will  generally  be  very  much  less  than  the  upper  bound  C  (e.g., 
the  graph  in  Figure  6  which  required  79  comparisons  as  opposed  to  several 
million  for  the  complete  symmetric  graph  (with  loops)  on  14  vertices). 
However,  since  the  non-zero  elements  of  P  will  generally  disappear 
quickly  from  P  as  k  increases,  the  number  of  basic  operations  required 
to  carry  out  Danielson's  algorithm  will  often  approach  D.   To  illustrate, 
we  will  compare  the  number  of  basic  operations  required  to  analyze  the 
graph  in  Figure  17(a)  using  Danielson's  algorithm  with  the  number  of 
basic  operations  required  to  analyze  the  graph  in  17(b)  using  our  algori- 
thm. 


Figure  17      (a)  An  Undirected  Graph   (b)  A  Digraph 

Using  our  algorithm  on  Figure  17(b),  we  obtain  the  attainability  tree 
shown  in  Figure  18  (the  only  one  for  this  graph,  since  all  vertices 
appear  in  it) . 
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r3 

Figure   18.      The  Attainability  Tree   for   the  Graph   of  Figure    17(b) 
We   see   that   this   graph   requires   22  comparisons   of   integers. 

Using  Danielson's   algorithm  on   the   graph   in  Figure   17(a),   we   have 


and 


p2    3 


P2   = 


T1**3     V1 

V1    V1*!2| 


P~    - 


2  3    3  2 
V1V2*V3^3V2 


Vi^Wl     V2V1+V3Vl+r3T2 


V371 


■ 1   2     13     2  3       13     2  12  3 


V1V2^3V1^3V2 
TlWl 


We  see  that  this  required  78  basic  operations. 
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VI.   CONCLUSIONS 

We  have  developed  an  algorithm  for  locating  the  simple  cycles  of  a 
digraph  which  appears  to  be  superior  to  previous  algorithms  in  terms  of 
the  number  of  basic  operations  required  and  which  is  conservative  in  its 
use  of  memory  space  on  the  digital  computer.   In  particular,  we  have  analyzed 
our  algorithm  and  the  one  originated  by  Danielson  [3],  and  found  that  the 
upper  bound  for  the  number  of  basic  operations  using  Danielson 's  algorithm 
was  given  by 


D-   SCkn^-l)  , 


k=2 

while  the  upper  bound  on  the  number  of  basic  operations  using  our  algori- 
thm was  given  by 

C 


=  Zjkn(n-l)!/(n-k)! 


k=l 

k-1 
Since  n    is  greater  than  (n-1) .'  /(n-k)  .' ,  we  found  that  C  was  less  than  D. 

Moreover,  since  the  basic  operations  required  by  Danielson 's  algorithm 

are  more  complex  than  those  required  by  our  algorithm,  we  see  that  our 

method  represents  a  considerable  improvement  over  the  matrix  technique. 

By  introducing  the  function  p  defined  in  Chapter  IV,  we  were  able 
to  modify  our  algorithm  for  application  to  coalesced  graphs.   Thus,  we 
can  analyze  extremely  cumbersome  graphs  using  the  techniques  developed, 
without  being  required  to  perform  an  undue  number  of  basic  operations. 

It  is  our  opinion  that  we  have  accomplished  the  goal  we  set  out  to 
achieve.   However,  it  is  believed  that  the  algorithm  might  be  made  even 
more  efficient.   It  was  noted  after  the  proof  of  Theorem  II  that  some 
simple  cycles  may  be  located  more  than  once.   We  suggest  that  there  may 
be  a  way  of  eliminating  this  ambiguity.   If  this  could  be  done, then  the 
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number  of  basic  operations  required  for  certain  graphs  would  be  cut  down 
considerably. 

It  was  noted  that  Danielson's  algorithm  was  designed  for  undirected 
graphs.   Clearly,  it  is  not  very  efficient.  We  suggest  that  a  method 
similar  to  the  one  we  have  developed  could  be  devised  for  application 
to  undirected  graphs  which  would  be  superior  to  Danielson's. 
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